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We discuss the recent progress made in understanding how some fundamental resuhs valid in the 
quenched infrared sector of non-Abelian Yang-Mills theories generalize to the unquenched case. 
In particular, we derive in the continuum an equation that allows to predict how the momentum 
dependence of the quenched gluon propagator is affected by the presence of a small number of 
dynamical quarks. In addition, we present the most recent lattice simulations of the (Landau 
gauge) gluon and ghost two-point sector obtained using gauge field configurations with two light 
and two light plus two heavy twisted-mass quark flavors. A comparison between the results 
obtained by these two complementary methods is then carried out. 
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Our understanding of the infrared (IR) sector of non-Abelian Yang-Mills theories has consider- 
ably improved over the last few years. In particular, ab-initio quenched simulations [|l|,^ performed 
on large volumes lattices have unequivocally demonstrated that the (Landau gauge) gluon propaga- 
tor and ghost dressing function saturate in the (deep) IR region. Within the continuum formulation 
of the theory, these lattice results are in agreement with, e.g., the solutions of the corresponding 
all-order Schwinger-Dyson equations (SDEs) [||, ||] and exact renormalization group equations [||]. 
This has caused a paradigmatic shift among practitioners: the gluon is now thought to acquire a 
momentum-dependent mass m{q^) whose magnitude can be large at IR momenta, but rapidly van- 
ishes with increasing spacelike momenta {i.e., » Aqj,^), thereby maintaining full accord with 
perturbative QCD. 

Here we report on our recent studies on the extension of these quenched results to the case in 
which dynamical quarks are present, both in the continuum as well as on the lattice. In particular: 

• From the continuum side, we describe, within the general formalism based on the pinch 
technique (PT) [g 0, H |, [T^ and the background field method (BFM) and known in 



the literature as the PT-BFM framework ^% pjj |T^, how one can predict the effect that 
the presence of a small number of dynamical quarks will have on the shape of the gluon 
propagator, starting from the quenched results [|5|] ; 

• From the lattice side, we will present the most recent data obtained in [|^] where the full 
non-perturbative ghost and gluon two-point sector was computed by using gauge field con- 
figurations (gauge fixed in the Landau gauge) with two light and two light plus two heavy 
twisted-mass quark flavors. 

A preliminary comparison between the results obtained will be finally presented. 



1. Schwinger-Dyson results 

In the PT-BFM framework it is possible to develop an approximate method for computing 
nonperturbatively the effects on the gluon propagator induced by the presence of a small number 
of dynamical quarks [p^]. 

To see how this can be done, let us start by observing that when dynamical quarks are included 
in the calculation of a particular Green's function, there are two effects that naturally occur: (/) a 
new set of diagrams involving quark loops explicitly appears in the unquenched version of the 
corresponding SDE; (//) the original quenched graphs get also modified, given that the presence 
of the quarks will in turn modify propagators, vertices and kernels that appear inside them. In 
particular, specializing to the case of the (PT-BFM) SDE describing the gluon propagator, one 
observes that in the full QCD case the original diagrams constituting the quenched version of the 
SDE, get supplemented with an extra graph containing the quark loop; however, the inclusion of 
this extra term will also affect indirectly the original (quenched) graphs since nested quark loops 
will now appear in them. 

To proceed further, one then assumes that the quark-loop diagram captures the leading cor- 
rection term, while additional corrections stemming from diagrams with nested quark loops are 
subleading. Neglecting these latter terms in front of the former, allows therefore to express the 
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unquenched propagator as the sum of two contributions: the quenched propagator we started from 
plus the quark loop diagram. Then, in the PT-BFM framework the unquenched gluon propagator 
Ag(^^) can be obtained from the quenched propagator A{q^) through the master formula (Euclidean 
space) |[T|] 

A.tf) = ^.^ ^ . (1-1) 



1 + 



Let us briefly recall what the different terms appearing on the right-hand side of the equation above 
correspond to. 

To begin with, X{q^) represents the PT-BFM scalar cofactor resulting from the calculation of 
the quark loop diagram; specifically, one has 

X{q^) = -^^ jTv[fS{k)f^{k + q-k-q)S{k + q)\ , (1.2) 

where the J-dimensional integral (in dimensional regularization) is defined according to 
= J df^k, with d = 4 — e and jJ. the 't Hooft mass. In addition, S represents the full quark 
propagator, defined as 

(p) = -i [A{p)^ - B{p)] = -iA{p) \i> - J{{p)\ . (1.3) 

The ratio ^{p) = B{p) /A{p) represents the dynamical quark mass function; finally, the vertex Fy 
is the full PT-BFM quark-gluon vertex, which, due to its characteristic Abelian-Uke nature, satisfies 
the linear Ward identity (WI) 

ip^^f^{p,,P2,p^)=S-\p,)-S-\-p2). (1.4) 

The appearance of vertices satisfying linear Wis as opposed to the usual nonlinear Slavnov-Taylor 
identities (STIs) satisfied by conventional vertices, represents the distinctive feature of the PT- 
BFM framework, and simplifies the problem considerably, as it allows to employ for such vertices 
Abelian-like Ansatze. 

Next, the function G{q^) is a special Green's function particular to the PT-BFM which achieves 



the conversion from the PT-BFM gluon two-point sector to the conventional one [13, 14|. In the 
Landau gauge G{q^) is known to coincide with the Kugo-Ojima function |17, 18]; in addition. 



this function can be traded for the ghost dressing function F thanks to the exact relation [ |17| , |18| ] 
F^\q'^) = 1 +G{q^)+L{q^), which, due to the fact that one has L{q^) < G{q^) in the entire 
momentum range [|l7|], lead to the approximate result F^^{q^) ?a \ + G{q^). 
Finally, the quantity X^{q^) is defined as 

X\q^)=ml{q^)-m\q^), (1.5) 

and corresponds to the difference between the dynamically generated gluon mass in the unquenched 
and the quenched cases, respectively. Evidently, the master equation ( |1.1D has been derived under 
the additional assumption that the inclusion of light quark flavors (say, e.g., an up and down type 
quarks, with constituent masses of about 300 MeV) will affect but not completely destabilize the 
mechanism responsible for the generation of a dynamical gluon mass, and that their effect may 
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be considered as a "perturbation" to the quenched case (the inclusion of loops containing heavier 
quarks has been shown to give rise to numerically suppressed contributions, consistent with the 
notion of decoupling [p^]). 

One can next use the quenched SU (3) results obtained in large volume lattice simulations [||] 
in order to evaluate the propagator A{q^) and the ghost dressing function F{q^) [and therefore 
indirectly the special function G{q^) due to the relation between them]. Thus, what is left to be 
computed are the quantities X{q^) and one proceeds then as follows: 

• As far as the quark loop X{q^) is concerned, its evaluation comprises two steps. First, the 
nonperturbative behavior of the functions A{p) and ^{p) appearing in the definition of the 
full quark propagator is obtained by solving the quark gap equation. Particular care is needed 
for approximating the non-Abelian quark-gluon vertex appearing in this equation, since 
the latter is the conventional vertex (and not the PT-BFM vertex) , therefore satisfying 
the usual STI involving the ghost dressing function as well as the quark ghost scattering 



kernel []19|]. It turns out that the inclusion of this dependence in the vertex Ansatz used 
for solving the resulting nonlinear system of integral equations, is crucial in order to boost 
up the dynamically generated quark masses to phenomenologically acceptable values [pO|]. 
Once the functions A{p) and ^{p) are known, one can proceed to the calculation of the 
scalar factor X{q^), employing for the PT-BFM quark-gluon vertex F^ a suitable Abelian- 



like Ansatz, e.g., the Ball-Chiu (BC) vertex [ ]21[ ] or the Curtis-Pennington (CP) [gj] vertex. 



Irrespectively from the vertex chosen, one has that in the q ^0 limit X (0) = [|15|] ; thus, the 
presence of fermions will not affect directly the saturation point of the (unquenched) gluon 
propagator (we will return on this important issue in the following point). 

As for the quantity X^{q^) its first principle determination requires the knowledge of the 
equation which governs the dynamics of the effective gluon mass. This equation was derived 



in [ ]23| , |24j] and turns out to be quite complex; in particular, it depends on the full gluon prop- 
agator: thus even though X(0) = 0, the unquenched saturation point will be different from 
the quenched one, since the overall shape of the propagator will be affected by the presence 
of dynamical quarks. The determination of the full ^^{q^) using the aforementioned equa- 
tion deserves a separate thorough study; in the ensuing analysis we will instead follow [p^], 
restricting ourselves to extracting an approximate range for A^, by employing a suitable 
extrapolation of the (unquenched) curves obtained from intermediate momenta towards the 
deep IR (thus we will effectively treating it as a constant). 

The results of the procedure so far described are shown in Fig. |l|, and can be summarized as 
follows. The basic effect of the quark loops (two families with a constituent mass of the order of 
300 MeV) is to suppress considerably the gluon propagator in the IR and intermediate momenta 
regions, while the ultraviolet tails increase, exactly as expected from the standard renormalization 
group analysis. In addition, the extrapolation procedure used (a cubic B-spline method applied 
starting from different momentum values) robustly predicts that the inclusion of light quarks makes 
the gluon propagator saturate at a lower point, which can be translated into having a larger running 
gluon mass. As far as the gluon dressing function q^A{q^) is concerned, one observes a suppression 
of the intermediate momentum region peak, a feature already reported in early lattice studies. 
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Figure 1 : (Top panels) Solution of the quark gap equation: A^^{p) (left) and dynamical quark mass 
■^{p) (right) renormalized at /i = 4.3 GeV; dotted black curves correspond to the BC vertex, while 
dashed blue curves to the CP vertex. {Bottom panels) Comparison between the quenched and the 
unquenched gluon propagator (left) and dressing function (right), obtained using the BC vertex. 
The shaded striped band in the left plot shows the possible values that ^q{Q) can assume at zero 
momentum depending on the extrapolation point used; in the case of the dressing function (which 
is basically insensible to the IR saturation point) we used a curve with an extrapolation point at 
= 0.05 GeV^. The quenched lattice results of [||] are also displayed for comparison. 



2. Lattice results 

From the point of view of lattice simulations, a comprehensive quantitative study of the gluon 
an ghost sectors in SU(3) QCD incorporating the effects stemming from the presence of two light 



and two light and two-heavy mass-twisted lattice flavors has been presented in [|16|] (see references 
therein for earlier lattice studies). The degenerate light quark masses used for the simulations 
ranged from 20 to 50 MeV, while the strange quark was roughly set to 95 MeV and the heavy 
charm to 1.51 GeV (in MS at a scale of 2 GeV). These values corresponds to a mass for the lightest 
pseudoscalar in the range approximately between 270 and 510 MeV; finally, the biggest volume 
simulated corresponded to an asymmetrical box of roughly 3^ x 6 [fm^]. 

The results obtained (for an arbitrary gauge copy selected through a gauge fixing algorithm 
using a combination of over-relaxation and Fourier acceleration) are shown in Fig. ^ and are in 
full agreement with the SDE study discussed in the previous section. Indeed, one sees again that 
the effect of the presence of dynamical quarks on the gluon propagator A is twofold: a suppression 
of both the "swelling" region at intermediate momenta and the saturation value in the deep IR. In 
addition, one observes that the size of the effect increases as more light flavors are added. 
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Figure 2: (rop panels) Landau gauge lattice results for the gluon propagator (left) and dressing 
function (right) for two light degenerate quarks (blue up triangles and green down triangles) and 
two light plus two heavy quarks (orange squares and red circles). (Bottom panels) Lattice results 
for the ghost dressing F{q^) and G{q^) functions (color code as above). In all cases the data have 
been renormalized at /i =4.3 GeV. 



On the other hand, the effect on the ghost dressing function F is much milder and is diametri- 
cally opposed to the one encountered for the gluon case, i.e., it consists in a small increase of the 
saturation point. Since in the SDE for the ghost, dynamical quarks enters only indirectly, either 
through the gluon propagator or via higher loop corrections to the gluon-ghost vertex, these results 
give support to the assumption (done in the SDE analysis of the previous section) that quark cor- 
rections due to nested quark loops can be expected to be small; in addition, it also shows that the 
approximation of the unquenched F{q-^) and G{q^) with the corresponding quenched functions is 
valid to a very high degree of accuracy. 

3. Comparison and outlook 

In Fig. |3| we plot a comparison between the lattice results and the SDE computations for the 
case of two degenerate light quarks. The figure shows an excellent qualitative agreement and a 
rather favorable quantitative agreement even at this level of approximation. 

The discrepancies observed are due to the way in which has been determined in [p^. 
On the one hand, the extrapolation procedure overestimates the saturation point of the quenched 
propagator (notice that the blue lattice data point at zero momentum is affected by large volume 



6 



Unquenching the infrared sector of QCD 



Daniele Binosi 



5 

^ 4 



> 

o 



Jl=4.3 [GeV] N(=2, F1=3.M, V=24,x48 ^ 
N(=2, P=4.2D, V=48 x96 ^ 
SDE, N,=2 . 



g2 [GeV^l 





; IGeVl 



Figure 3: Comparison of SDE and lattice results for the case of two light fermions. The extrapo- 
lation of the SDE curve towards the IR is obtained starting from = 0.05 GeV^. 



artifacts, see [16]); on the other hand, the lack of knowledge of the full momentum dependence 
of implies that the SDE curve (as is best seen in the plot of the dressing function) does not 
reproduce quantitatively the data in the intermediate momentum region (since the dynamical gluon 



mass is dropping pretty fast with [24], this is precisely the region in which we expect the biggest 
difference between the SDE and the lattice results). The determination of the full from the 



recently derived all-order mass equation [24] represents clearly the next step for this approach. 

From the lattice point of view, on top of the obvious improvements which can be obtained by 
running simulations with bigger volumes and statistics and smaller pseudoscalar masses, it would 
be certainly very interesting to implement directly the simulations of PT-BFM Green's functions 
using the canonical transformation technique recently developed in 
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